Hartree-Fock-Bogoliubov Theory of Polarized Fermi Systems 
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Condensed Fermi systems with an odd number of particles can be described by means of polar- 
izing external fields having a time-odd character. We illustrate how this works for Fermi gases and 
atomic nuclei treated by density functional theory or Hartree-Fock-Bogoliubov (HFB) theory. We 
discuss the method based on introducing two chemical potentials for different superfluid compo- 
nents, whereby one may change the particle-number parity of the underlying quasiparticle vacuum. 
Formally, this method is a variant of non-collective cranking, and the procedure is equivalent to the 
so-called blocking. We present and exemplify relations between the two-chemical-potential method 
and the cranking approximation for Fermi gases and nuclei. 

PACS numbers: 03.75.Mn, 05.30.Fk, 21.60.Jz, 31.15.cj 



I. INTRODUCTION 

Due to the experimental realization of strongly in- 
teracting atomic Fermi gases, there has been increased 
theoretical interest in properties of asymmetric two- 
component superfluid Fermi systems with unusual pair- 
ing configurations 0, 0, i i, 0, 1, ©, M EH- Of 
current interest are the properties of spin-polarized con- 
densates having an unequal number of spin-up and spin- 
down fermions. One of the condensation possibilities is 
the "breached-pair" (BP) superfluid state, in which the 
normal phase coexists with a full pairing of the minor- 
ity fermions at high spin imbalances. Recently, a gap 
in a single-particle excitation spectrum of a highly spin- 
imbalanced sample has been observed experimentally in 
a 6 Li atomic condensate fl2^ . 

Atomic nuclei can also exhibit interesting pairing prop- 
erties, including BP superfluidity, although the fraction 
of the polarization is typically quite small. Examples 
include: 

• Pairing isomers [HI, EH, EH having an appreciably 
smaller value of the pairing gap than the ground 
state. In the roots of pairing isomerism is very dif- 
ferent single-particle level density in the vicinity of 
the Fermi surface in coexisting configurations, usu- 
ally associated with different nuclear shapes. Such 
states have been found in, e.g., ^4«80, 100, and 190 
mass regions. 

• Gapless superconductivity at high spins of atomic 
nuclei [H, Here, quasiparticle energies in the 
rotating frame become negative at certain rota- 
tional frequencies; i.e., there is no gap in the quasi- 
particle spectrum. However, the pairing energy, 
or the expectation value of the pairing condensate, 



still remains large. 

• Angular-momentum aligned nuclear configurations, 
such as high-spin isomers, in which pairing cor- 
relations ma y b e reduced due to blocking (see 
Refs. [H,[H)[13l and references quoted therein). 

In this paper we will briefly rederive some of the 
generic results of the Hartree-Fock-Bogoliubov theory 
(also called the Bogoliubov de-Gennes theory) for po- 
larized quasiparticle states, using constrai ning fields to 
reach the states, as was done in Refs. @, flO| to study 
the polarized atomic condensates. We shall call this the 
"two- Fermi level approach" (2FLA). In particular, odd- 
mass systems require special attention in the HFB theory. 
In the usual blocking approach, the selected quasiparti- 
cle state becomes occupied and this requires a modifi- 
cation of the HFB density matrix. In 2FLA, the low- 
est quasiparticle orbit becomes occupied by changing the 
particle-number parity of the vacuum through the exter- 
nal field without modifying any of the vectors explicitly. 
We show that for spin systems (e.g., Fermi gases) the 
2FLA is equivalent to the standard rotational cranking 
approximation while this is not the case for atomic nuclei, 
in which polarization is due to the angular momentum 
alignment. In both cases, however, 2FLA can be viewed 
as a "vacuum selector" by means of the non-collective 
cranking. 

The paper is organized as follows. Section |TT] briefly 
discusses the concept of Bogoliubov quasiparticles. Par- 
ticular attention is paid to the choice of the HFB vacuum, 
the way the particle-number parity is encoded in the Bo- 
goliubov matrix transformation, and self-consistent sig- 
nature symmetry of HFB and its relation to time rever- 
sal. The HFB extension to the case of two-component 
systems (2FLA) is outlined in Sec. IIII1 and its relation 
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to non-collective cranking and the blocking procedure is 
discussed. Section [TV] shows numerical examples, both 
from atomic and nuclear physics, based on the HFB ap- 
proach. Finally, Sec. fVl contains the main conclusions of 
this work. 



The HFB vacuum |\E') is a zero-quasiparticle state: 



=0. 



(3) 



II. THE QUASIPARTICLE FORMALISM 

We begin by recalling basic equations of the quasiparti- 
cle formalism, which historically is attributed to Gor'kov, 
Bogoliubov, and de Gennes. While these equations and 
definitions are admittedly very well known, there are sev- 
eral aspects of the quasiparticle approach that are seldom 
discussed; hence, they are worth bringing to the atten- 
tion of a wider community. We shall discuss these lesser- 
known aspects in the subsections following the general 
introduction to HFB. 

The HFB wave functions are quasiparticle product 
states. The quasiparticle annihilation operators a M are 
defined as linear combinations of particle annihilation 
and creation operators by the Bogoliubov transforma- 
tion, 



(1) 



The matrices U and V satisfy the following canonical 
conditions: 



Complete information about \^>) is, in fact, contained in 
the generalized density matrix 1Z, 



PuH , 

K + 



which in terms of matrices U and V reads: 



K = 



*T T T 
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(4) 



(5) 



The variational principle implies that the self-consistent 
density matrix TZ commutes with the quasiparticle 
Hamiltonian 



H 



h-XI 
A+ 



A 

-h* + XI 



(6) 
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(2a) 
(2b) 
(2c) 
(2d) 



where h = t + Y is the single-particle Hamiltonian and T 
and A are particle- hole and particle-particle mean-fields, 
respectively. The HFB equations can be written in a 
matrix form: 



h-XI , A 
A+ , -h* + XI 



U , V* 
V , U* 



r 



u 
v 



V* 

u* 



E 




-E 



(7) 



where E is a diagonal matrix of quasiparticle energies 
Ep. Columns of eigenvectors, 



y* 

u* 



x ■-- 



(8) 



are called occupied and empty quasiparticle states, re- 
spectively, because they are eigenvectors of the projective 
matrix TZ with eigenvalues 1 and 0. 



A. Choice of occupied quasiparticle states 

In practical applications, one often assumes that the 
quasiparticle energies in matrix E are all positive. How- 



ever, this is not at all required by the variational princi- 
ple. Indeed, while eigenvalues of H must come in pairs 
of opposite quasiparticle energies (E^, —E^), the theory 
says nothing on whether a positive or a negative one en- 
ters matrix E. 

This question is evidently related to the mundane 
problem of the order in which we arrange eigenvectors 
of H in Eq. ((7|). Any such order is allowed, provided 
one eigenvector of the pair (E^, —E^) is put into the first 
half of the spectrum (empty quasiparticle states), and 
the other one is put in the second half of the spectrum 
(occupied quasiparticle states). However, which one goes 
where is in principle arbitrary. The signs of quasiparti- 
cle energies in E are arbitrary too, and thus variational 
equations have not one, but many solutions. Naturally, 
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in practical applications, the choice of occupations is mo- 
tivated by physics as the energy of the system may be 
different depending on how the arbitrariness of choos- 
ing signs of quasiparticle energies is resolved. It is then 
plausible, and often assumed, that the lowest total energy 
is obtained by occupying states corresponding to nega- 
tive quasiparticle energies, and leaving empty those cor- 
responding to positive quasiparticle energies, i.e., by as- 
suming that positive quasiparticle energies are collected 
in matrix E. One should stress at this point, however, 
that there is no a priori reason why such a choice must 
guarantee obtaining the lowest total energy, and exam- 
ples to the contrary are available in numerous applica- 
tions. Thus, which quasiparticle states should be occu- 
pied is a matter of a specific physical situation, and not 
a rule cast in stone. 

The problem of choosing occupied quasiparticle 
states is particularly conspicuous when iterative (self- 
consistent) methods are used for solving the non-linear 
eigenproblem ||7J). The self-consistent procedure can be 
described as the sequence of the following steps: 

1. Find eigenstates of the quasiparticle Hamiltonian 
©■ 

2. Choose the occupied quasiparticle states. 

3. Calculate the generalized density matrix ([5]), i.e., 
the density matrix (p) and pairing tensor (k). 

4. Calculate the particle-hole (T) and particle-particle 
(A) mean fields to determine the Hamiltonian ([6]) 
in the next iteration. 

5. Return to step 1. 

Again, this self-consistent sequence of steps is described 
in all textbooks; however, the crucial step 2 is seldom ever 
mentioned. However, in many applications, this step is 
absolutely essential for obtaining a convergent algorithm. 
Moreover, choosing negative occupied quasiparticle ener- 
gies at each iteration may sometimes lead to divergent 
iterations. 

The problem here is in finding ways of identifying 
quasiparticle states (tags), which are independent of their 
energies, and which would allow for pinning down the 
structure of each individual state. Ideally, tags could be 
provided by quantum numbers of conserved symmetries. 
For example, spherical symmetry allows each state to be 
characterized by the standard quantum numbers nljm. 
Then, decisions of occupying quasiparticle states can be 
made based on lists of quantum numbers. 

A particular version of using conserved quantum num- 
bers as tags of quasiparticle states was pioneered in nu- 
clear physics under the name of non-collective rotation, 
or non-collective cranking approximation [III Gil HH > an d 
this method is a central theme of our discussion. It con- 
sists of using in Eq. ^ the single-particle Routhian h' 
instead of the single-particle Hamiltonian h, where 



and S is the operator of a conserved symmetry. This may 
look like an attempt of introducing a Lagrange multiplier 
to constraining the average value of S, but in fact it is 
not. Indeed, since quasiparticle states are eigenstates of 
S, the quasiparticle energies are simply reordered as a 
function of As without any modification of the structure 
of quasiparticle states or total energy. Using Routhian h! 
may thus facilitate selecting proper quasiparticle states 
into the set of occupied ones, depending on the physical 
situation corresponding to the chosen symmetry S. In 
short, non-collective cranking can be viewed as a config- 
uration selector. 

The method fails when quasiparticle states arc not 
eigenstates of symmetry operators, which can happen in 
unrestricted calculations. The practical solution consists 
of calculating overlaps of quasiparticle states with a fixed 
set of wave-functions and establishing in this way tags 
that are related to the closest similarities of structure 
between both sets. 

Only when a set of properly tagged quasiparticle states 
is occupied at step 2 of each iteration can one have a good 
chance of obtaining a converged solution. At the end of 
the day, one can check whether such a solution corre- 
sponds to the negative quasiparticle states being occu- 
pied or not. 

B. States with even and odd number of fermions 

The ambiguity of choosing the occupied quasiparticle 
orbits is particularly important for odd-particle systems. 
Equations H])-© look entirely the same irrespective of 
whether the state |^) represents even- or odd-particle 
systems. While the particle number N is not conserved 
by the product state \^f), i.e., this state is a linear su- 
perposition of components having different particle num- 
bers, the parity of the particle number, 7Tjv = (— 1) , is 
conserved. That is, the decomposition of |\&) may con- 
tain cither even- or odd-particle number components, but 
never both. 

Of course, information on whether the particular vac- 
uum |\&) is 7r/v-even or odd is contained in matrices U 
and V of the Bogoliubov transformation {TJ), i.e., it must 
depend on choices made for the occupied quasiparticle 
states. Specific choices of occupied quasiparticle states 
must be made in order to obtain even (ttn — 1) or odd 
(ttjv = — 1) states. 

In which way is the information on 7r/v encoded in 
matrices U and VI In the simplest situation, the set 
of occupied quasiparticle states are such that U is non- 
singular, det(U) =/= 0. Then we can determine the matrix 
Z = VU~ X , and express the vacuum \*&) through the 
Thouless theorem [17j : 

I*) = |$}cvc„ = AAexp f-i^Z+a+a+) I )' ( 10 ) 



ti := h - \ S S, 



(9) 



where the normalization factor is Af = det 1 ^ 2 (U + ). From 
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this expression we see that |\&) has even number parity, 

7TJV = 1. 

Let us now make another choice of selecting the occu- 
pied states: we replace one (/ith) column of the matrix 
of occupied states tp by the same (//th) column of the 
matrix of empty states x ( see Eq. (JH)), i.e., 

Uit = [ZK V f V + ^ (11a) 
I V u-u for v = fi, K ' 

V M _ / Vv'v for v + (uh \ 
v ' v ~ 1 Efr„ forr/ = M , lllDJ 

In this way, we replace the quasiparticle energy of the 
occupied state by — E^. We do not imply here that 
E^ must have been positive, so we have replaced a pos- 
itive quasiparticle energy by a negative one - we could 
have just replaced the negative one by the positive one. 
In fact, as discussed earlier, there is no such rule that oc- 
cupations leading to the original non-singular matrix U 
must correspond to all quasiparticle energies being posi- 
tive. 

Now it is a matter of simple algebra to see that for 
such a choice of occupied quasiparticle states, the vacuum 
state reads 

\*)&=Ma+exp (-hJ2 Z »y\ l°>> ( 12 ) 

\ flV J 

and is a manifestly odd state, ttn = — 1 5 which we call a 
one-quasiparticle excitation of |$) C vcn- 

Note that in the odd state (p~2|) , the matrix Z and nor- 
malization constant TV are defined through the original 
matrices U and V of the even state (|10p , and not through 
those after the column replacement as in Eq. (fTTj) . In- 
deed, it is easy to see that the Thouless theorem does 
not work for the one-quasiparticle states because matri- 
ces u\jfl are singular. This is obvious from Eq. (|2d[) . 
whereby each column of matrix V* is a linear combina- 
tion of columns of matrix U, i.e., 

V* = -U(V+(U T )- 1 ). (13) 

Therefore, after the column replacement as in Eq. (fTTj) . 

matrices of one-quasiparticle states, Uj^, become sin- 
gular and have null spaces of dimensions D—l. Con- 
sequently, the corresponding matrices 1 — p^' — 
U(p)jjM+ : cf Eqs . Q an d ©, have exactly one eigen- 
value equal to zero. Hence, the occupation numbers 
(eigenvalues of the density matrices pW) of one of the 
single-particle states are in each case equal to 1. This fact 
is at the origin of the name "blocked states" attributed to 
one-quasiparticle states (fT2")l . These states contain fully 
occupied single-particle states that do not contribute to 
pairing correlations. 

We can continue by building two-quasiparticle states 

|*>&£ = Na+a+ cxp ( -| ]T Z+,o+a+\ |0), (14) 



which are manifestly particle- number parity even, 7Tjv = 
1. They correspond to two columns in ip replaced by two 
columns of x, an d the corresponding matrices [/(w ) have 
null spaces of dimensionality D~0 or D=2. We now have 
a very definite prescription for telling which Bogoliubov 
transformations correspond to even and which to odd 
states, i.e., itn = ( — 1) D , where D is the dimensionality 
of the null space of matrix U. 

The main conclusion of this section is that vacuum 
states of given ttn are obtained by making appropriate 
choices of occupied quasiparticle states. In particular, 
one should begin by selecting one even state (fT0|) repre- 
sented by a non-singular matrix U, which one can call a 
reference state, and than proceed by building on it one-, 
two-, or many-quasiparticle excitations. Note that con- 
structing odd states is best realized by first building an 
even reference state and then making one-quasiparticle 
excitations thereof. This is best done by blocking spe- 
cific quasiparticles, i.e., replacing columns of matrices as 
in Eq. (|TT|) . After the self-consistent procedure is con- 
verged for each blocked state, one may select the lowest 
one as the ground state of an odd system, and consider 
the higher ones as good approximations of the excited 
odd states. It is obvious that self-consistent polarization 
effects exerted by blocked states, which will be taken into 
account by iterating Eq. l[7]). may render reference states 
of every blocked configuration to be different from one 
another. 

Note also that in the above analysis we did not talk 
about the average particle numbers, which can be even, 
odd, or fractional, depending on the value of the Fermi 
energy A in Eq. ([5]). Thus one can, in principle, consider 
odd states with even average particle numbers, or even 
states with odd average particle numbers. The latter ones 
provide especially useful reference states for building one- 
particle excitations on top of them, because they require 
the smallest readjustment of the Fermi energy between 
the reference state and a one-quasiparticle excitation. 



C. Signature symmetry 

In this section we apply methods of occupying quasi- 
particle states, outlined in Sec. Ill Al to a physical situa- 
tion where the system has a conserved signature symme- 
try, which can be used within the cranking approximation 
[l7l . [IH, [2l[ . The signature operation [22] is a rotation by 
7r around one direction in space, which is conventionally 
called the cc-axis: 

R x = exp(-?7rJ 2; ), (15) 

where J x denotes the total angular-momentum operator 
along the x axis. The signature operator is manifestly 
unitary, R X R X = 1. Since a rotation by 2ir reverses the 
phase of fermion wave functions, the square of the signa- 
ture operator gives the particle- number parity, R x = ttn ■ 
Therefore, R x is hermitian and antihermitian in even and 
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odd spaces, respectively. In particular, in the single- 
particle space, the signature is a unitary antihermitian 
operator. 

Since the signature R x and time-reversal T operators 
commute, 



T + R X T = R x , 



(16) 



signature is a time-even operator. Therefore, for non- 
rotating systems (i.e., without time-odd fields), signature 
is equivalent to the time-reversal symmetry T, but it is 
more convenient to use, because R x is a linear and not an 
antilinear operator 23]. For states with nonzero angular 
momentum, where T is internally broken, R x is often 
still preserved [23|, [24], HH . The single-particle (and one- 
quasiparticle) states may then be classified according to 
the signature exponent quantum number a (26| : 



R x \ak) 



\ak), 



(17) 



where a takes the values of ± 1 /2- For conserved signa- 
ture, the HFB mean fields h and A commute and anti- 
commute with R x , respectively, and the HFB equations 
([7]) can be written in a good signature basis: 

H a \afi) = E ai t\ct(i), (18) 

where the HFB Hamiltonian matrix ((SJ) in one signature 
reads 



k. — ak' 



j { _ | \h — X) a k : ak' 



(19) 



and the two-component quasiparticle wave function is 

U afl 



|a/x) = 



yap. 



(20) 



A quasiparticle state with good signature is a linear 
combination of states in time-reversed orbits [23| : 

|l/2,fc> = ^[-|fcr> fe )+ 7 r(-l)^- 1 /2T|fcr! fe )](21) 
| — 1/2, Af> = ^[f|Wi fc >+7r(-l) n *- 1 / 2 |fen fc >],(22) 

where is the eigenvalue of the z component of 
the single-particle angular momentum, J z , and we 
have adopted the phase convention according to which 

f\njfl) = TT(-iy +n \TTj,-n) 0, El, Where 7T = ±1 is 

the parity quantum number. If the Kramers degeneracy 
is present, the description in terms of Kramers doublets 
and signature doublets is equivalent. It is for polarized 
systems having time-odd mean fields that the use of the 
signature symmetry is superior. 

The HFB equation (fT8| has the (quasiparticle- 
quasihole) symmetry ([7]): for each state |a/z) of a given 
signature, there exists a conjugate state of opposite sig- 
nature | — afl), opposite energy: 



(23) 



and the quasiparticle wave function given by 

| Ck-jl) [ jjotp* 



(24) 



By this symmetry, one needs to solve the HFB equation 
(|18p only for one signature, obtaining positive and neg- 



ative quasiparticle energies E afl . Therefore, the entire 
set of negative quasiparticle energies is composed of two 
groups: (i) the negative ones E afl obtained directly from 
the HFB equation solved for signature a, and (ii) the in- 
verted positive ones (f23|) . which correspond to states of 
signature —a. 

The zero-quasiparticle HFB reference state (fTU)) , rep- 
resenting the lowest configuration for a system with even 
number of fermions, corresponds to a filled sea of Bogoli- 
ubov quasiparticles with negative energies (Fig.^a)). In 
a one-quasiparticle state, representing a state in an odd 
nucleus, a positive-energy state is occupied and its con- 
jugated partner is empty (Fig. [T{b)). 



E- 2 \ 



-Ej 
-E, 



(a) 



0-qp (vacuum) 
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c 2 

E] 
E] 


-Er 
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-O- 



1-qp 

a=1/2 



Ei 



-Ej 
-E T 
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-o- 

-O- 



2-qp 

a=0 



FIG. 1: (Color online) Quasiparticle content of three config- 
urations: (a) vacuum; (b) the lowest one-quasi-particle state 
with a=l/2, accessible via 2FLA; (c) the lowest two-quasi- 
particle state with a=0, not accessible via 2FLA. 

The exchange of the eigenvectors (U^, V^) and 
(V^* ,C/ M *), that have opposite signatures, corresponds 
to the exchange of columns in the <p and x matrices dis- 
cussed in Sec. Ill Al and reverses the particle-number par- 
ity 7Tjv- The density matrix and a pairing tensor of a one- 
quasiparticle state (fT2"|) can be obtained from Eqs. (TTTj) 



Potk,a>l 
Z ak,-al 



f ak,al 
K ak,-al 



V -al V -ak 



u -al u -ak' 



(25a) 



where p° and n° correspond to the reference state (|10[) . 



III. TWO FERMI LEVEL APPROACH 

The main idea behind the 2FLA 0, [l(| is to force a 
nonzero spin polarization in the system by finding the 
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ground-state condensate in the presence of an external 
field that favors one spin over the other. In the lan- 
guage of signature, "spin-up" corresponds to a = 1/2 
while "spin-down" corresponds to a = —1/2. The po- 
larization is achieved by adding to the Hamiltonian the 
single-particle field 

h s = i\ s R x , (26) 

i.e., constructing the Routhian ([9]), h! — h — h s . The 
imaginary unit must be put in the definition of h s , be- 
cause in the single-particle space R x is antihermitian. 
Consequently, h s is time-odd, cf. Eq. ([TO)) . Since R x and 
the quasiparticlc Hamiltonian commute, adding h s repre- 
sents a non-collective cranking; hence, the quasiparticle 
routhians must be linear in X s . 

The field (T2"6"l) will raise the Fermi energy A1/2 of the 
subsystem having signature a = +1/2 by an amount 
X s and lower the Fermi energy A_i/ 2 of the a = —1/2 
subsystem by the same amount. The relations between 
chemical potentials in 2FLA read: 

X a =\ + 2aX s , (27) 

where 

^ = 2 (^1/2 + A-1/2) , A s = 5 (A1/2 — A_i/2j . (28) 
The HFB Routhian matrix of 2FLA can be written as 



H S =H- 2a\ s 2, (29) 

where H is the matrix (fT9|) corresponding to X s = and 
X is the unit matrix. 

Since the added term is proportional to the unit ma- 
trix, its only effect is to shift the HFB eigenvalues up or 
down, 

K» = E a ^-2a\ s , (30) 

where E afl are the eigenvalues of TL. Therefore, when 
plotted as a function of A s , the energies are straight 
lines with slopes 2a = ±1, as schematically depicted in 
Fig. [21 The Hamiltonian (TT9"]) usually represents an un- 
polarized system, in which case the quasiparticle energies 
E afl are degenerate (Fig. [2J top). If TL has time-odd fields 
(due, e.g., to an external magnetic field or nonzero an- 
gular velocity), this Kramers degeneracy is lifted (Fig. [21 
bottom). 

A. Interpretation of polarizing field h s 

The external field (|2"6"1) has a particularly simple inter- 
pretation for spin systems that are spherically symmetric 
in space (e.g. in a spherically-symmetric trap). Since in 
this case we can neglect the orbital part L of the angu- 
lar momentum in J = L + s, the polarizing field can be 
written as 

h s = i\ s exp(-iTrs x ) = 2X s s x , (31) 




FIG. 2: (Color online) One-quasiparticle levels of both signa- 
tures (a— 1/2: solid line; a=-l/2: dotted line) as functions 
of A s for a HFB Hamiltonian (|19l) without (top) and with 
(bottom) time-odd fields. In the latter case, the Kramers de- 
generacy at A s =0 is lifted, i.e., states \afi) and T\ajl) have 
different energies. The negative energy levels occupied in the 
quasiparticle vacuum are marked by thick lines. At the points 
A 3 ™' marked by stars, the quasiparticle level E n with a=l/2 
becomes occupied. If this level is not degenerate, or its de- 
generacy is an odd integer, the particle-number parity ttn of 
the vacuum changes as indicated. 



and 2FLA is equivalent to the (one-dimensional) rota- 
tional cranking approach [Tt], [3 ; in which the total en- 
ergy of the system at a fixed value of the total angu- 
lar momentum J x — (J x ) is obtained by adding to the 
Hamiltonian the one-body field —luJ x (rotational crank- 
ing term). Therefore, for spin systems, the angular ve- 
locity lu is simply equal to 2A S . 

A different situation occurs in atomic nuclei, where 
the spin degeneracy is lifted by the spin-orbit interac- 
tion. Then, projection of the single-particle angular mo- 
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mentum VL X is a good quantum number. Consequently, 
the eigenvalues of the signature operator (fT5)) are sim- 
ply e-"^-, i.e., a=l/2 for = 1/2, -3/2, 5/2, •• ■ and 
a=-l/2 for fi x = -1/2,3/2, -5/2, • • •. In this situation, 
there is no difference between signature and spin pro- 
jection, and odd states can indeed be obtained within 
the 2FLA. However, in this case, the angular momen- 
tum polarization (or spin polarization in a deformed trap 
rotating along the symmetry axis) is best modeled by 
the cranking approximation, where J x , rather than sig- 
nature R x , is used to build Routhian This is so, 
because for systems such as nuclei that are governed by 
j — j coupling, the polarizing field 2X s a does not distin- 
guish between individual single-particle alignments, i.e., 
between states with large fl x values that predominantly 
contribute to the angular momentum alignment and low 
Q x states that weakly respond to rotation. In this re- 
spect, the standard non-collective cranking approach has 
a certain advantage. 

Within the non-collective cranking approach, the 
single-particle cranking term becomes — u>Q x and the 
single-quasiparticle energies are linear in lo with slopes 
distinguishing between angular momentum projections. 
The corresponding quasiparticle vacua represent "opti- 
mal" states with maximally aligned angular momentum 
(maximum polarization). The minimization of the to- 
tal energy for intermediate values of angular momentum 
can be done by considering the lowest particle-hole ex- 
citations across the Fermi surface. These states are not 
accessible within the 2FLA, and they must be obtained 
by explicitly blocking quasiparticle states. 



B. Particle-number parity of the HFB vacuum 

Let us first analyze the 2FLA in the case when lev- 
els depicted in Fig. [5] are not degenerate. At low val- 
ues of A s , the quasiparticle vacuum corresponds to a 
system with even number of fermions. It is seen that 
at \ s =E\/2 t i the down-sloping lowest quasiparticle with 
a=l/2 crosses zero and becomes negative. Beyond that 
point, the HFB vacuum has one quasiparticle state occu- 
pied, as in the middle panel of Fig. [TJ Here, the particle 
number parity ttn changes from even (+1) to odd (-1), 
as discussed in Sec. Ill B | a nd also in the context of nu- 
clear rotations in Refs. [17l. l3l| . This can also be derived 
in a straightforward and explicit way by calculating the 
expectation value of the number-parity operator e l7r ™ in 
the HFB vacuum. For each subspace in the canonical 
representation, the number operator can be expressed 
7Tjv a =e i7r ( ftl +™ 2 ) = (1 - 2ni)(l - 2n 2 ), where 1,2 label 
the quasiparticle transformations for the negative eigen- 
values. The ground-state expectation values of hi, hi, 
and hih 2 are then evaluated in the usual way by ex- 
panding in the quasiparticle basis, normal ordering, and 
extracting the zero-quasiparticle term. 

By using the modified density matrices in the HFB 



equations and in the particle number equation for A, 

N = Tr(p), (32) 

one formally recovers the standard blocked HFB equa- 
tions for the lowest 1-qp (a=l/2) state. At still higher 
values of A s , the particle- number parity changes again 
at X s =E 1 / 2 ,2 when the second lowest quasiparticle with 
a=l/2 crosses zero. The associated two-quasiparticle 
configuration in a nucleus with N = Nq + 2 has the sig- 
nature index a=l. It is immediately seen that the lowest 
two-quasiparticle a=0 configuration of Fig. QJc), associ- 
ated with the so-called S-band in rotating nuclei, cannot 
be reached within the standard 2FLA. 

Following the discussion in Sec. IIII A[ it is worth not- 
ing that the choice of angular momentum quantization 
implies a different character of the angular alignment as- 
sociated with the change in the quasiparticle vacuum. 
In the case of ^-quantization and the absence of time- 
odd fields in the HFB Hamiltonian, 2FLA treatment of 
systems with odd particle number is equivalent to the so- 
called uniform filling approximation, in which a blocked 
nucleon is put with equal pro bability in each of the degen- 
erate magnetic substates |32l . |33| • It is only in the regime 
of non-collective rotation in which the angular momen- 
tum is quantized along the x-axis that the dynamics of 
angular momentum alignment can be properly treated 
and the full alignment can be reached. 

In cases when in the unpolarized system the Kramers 
degeneracy is the only one, the parity 7r/v of states ob- 
tained by occupying the lowest < quasiparticles 
do change at points where E^^ — 0, as indicated in 
Fig. [2] However, if apart from the Kramers degeneracy 
there is an additional two-fold, four-fold, etc. degeneracy 
of quasiparticle levels, all such states will be 7T7v-even. 
Therefore, in such situations, the 2FLA would fail to pro- 
duce odd- A systems as HFB ground states. 

For spherical spin systems, each level in Fig. [2] is 
(2£+l)-degenerate, which is an odd number, and the 
number parity does change at points = 0. However, 
the corresponding HFB vacuum represents a (2^+1)- 
quasiparticle excitation and not the one-quasiparticle ex- 
citation of Eq. (fT2"l) . This is so even if the average Fermi 
energy A is adjusted to have (on average) only one par- 
ticle more than that of the unpolarized system 0, EH • 
Such a situation corresponds to the so-called filling ap- 
proximation of orbitally degenerate quasiparticle states. 
Needless to say, the orbital filling approximation com- 
pletely neglects possible space-polarization effects that 
must, in principle, occur for true one-quasiparticle states. 

If the spin system has an axial symmetry in space (e.g., 
it is in an external axially symmetric trap), projection 
of the angular momentum on the symmetry axis A is a 
good quantum number. Moreover, the single-particle and 
quasiparticle states are then degenerate with respect to 
the sign of A. Here, each level in Fig. [2l except for A=0, 
is doubly degenerate. Therefore, none of the A>0 states 
obtained by occupying the < quasiparticles has 
7Tjv= _ 1, that is, odd particle number. In this case, states 
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with 7Tat=-1 cannot be obtained within 2FLA, and ex- 
plicit treatment within the blocking approximation, de- 
scribed in Sec. Ill CI is mandatory. 



IV. ATOMIC AND NUCLEAR HFB 
CALCULATIONS 

The first numerical example of 2FLA deals with a two- 
component polarized atomic condensate in a deformed 
harmonic t rap using the superfluid local density approxi- 
mation [l(| HI]. The system is described by a local energy 
density 

r(r) 3(3^ 2 ) 2 / 3 p 5/3 (r) , „ , , l2 

£(r) = a u -Y+p u -± iQ P +g eff (r)\ K (r)\ 2 , 

(33) 

where the local densities p(r) (particle density), r(r) 
(kinetic energy density), and n(r) (pairing tensor) are 
constructed from the quasiparticlc HFB wave functions. 
The parameters a u , f3 u , and the effective pairing strength 
9eff( r ) nave been taken according to Ref. [10| . We as- 
sume that the external trapping potential can be de- 
scribed by an axially deformed harmonic oscillator with 
frequencies [ItJ 

w 2 ± (S) = u 2 (S) (l + , ^(6) = ul{5) (l ^ , 

(34) 

with 

uj (5) = Cj ^1 + ~<J 2 J . (35) 

As in Ref. [ici[, we put fio>o=l' 

The calculations were carried out for systems with 
iV=30 and 31 fermions in a spherical (<5=0) and deformed 
((5=0.2) trap. The HFB equations were solved by us- 
ing the recently developed axial DFT solver HFB- AX [35| . 
The results are displayed in Fig. 3. 

In the spherical case, Fig. 3 (top), the spin degener- 
acy is lifted by the polarizing field h s . However, as dis- 
cussed in Sec. lIIlBl the orbital (2f+l)-fold degeneracy is 
present. Consequently, after the crossing point, the vac- 
uum becomes a (2£i+l)-quasiparticle state, where l\ is 
the orbital angular momentum of the lowest quasiparticle 
level. For 7V=30, the lowest quasiparticle excitation is a 
p state, i.e., above the crossing point the local HFB vac- 
uum becomes a three-quaspiarticle state. At the crossing 
point, the self-consistent mean-field changes abruptly. In 
particular, the chemical potential moves up as the num- 
ber of particles increases by one, and the pairing gap 
decreases due to blocking. This produces a sharp discon- 
tinuity around the crossing point, which can be seen in 
all three cases presented in Fig. 3. 

The middle portion of Fig. 3 illustrates the deformed 
case. Here, quasiparticle states are labeled using the an- 
gular momentum projection quantum number A onto the 



N=30 N=31 




K 

FIG. 3: (Color online) Quasiparticle levels (q=1/2: solid 
line; a=— 1/2: dotted line) as functions of X 3 for a two- 
component polarized atomic condensate in a deformed har- 
monic trap using the superfluid local density approximation 
[ioL l34j |. The calculations were carried out for jV=30 (left) 
and N—31 (right). The top panels correspond to the spheri- 
cal case (5=0). Here the quasiparticle levels are labeled with 
the usual spectroscopic designation of orbital angular momen- 
tum i. Each line represents a set of (2^+l)-fold degenerate 
states. In the deformed case (5=0.2; middle panels) the lev- 
els with the opposite values of A>0 are two-fold degenerate. 
The corresponding values of |A| ere shown in the inset. The 
Kramers degeneracy can be lifted by adding a cranking term, 
—u>aIx to the Hamiltonian (bottom panels, — oja=0.05). Here, 
every level is labeled by the A quantum number, see the inset. 
In order to make the plot less busy, the quasiparticle levels 
originating from very excited spherical g, d, and s shells are 
not shown in the two lower panels. See text for details. 



symmetry axis of the trapping potential (rc-axis). Be- 
cause of the Kramers degeneracy, levels with ±A are de- 
generate. That is, except for A=0, each quasiparticle 
state is two-fold degenerate. In the case presented in 
Fig. 3, the two lowest levels have A=0, hence they are 
associated with one-quasiparticle excitations. The third 
state has |A|=1, and its crossing does not change the 
particle-number parity of the HFB vacuum. 

The two-fold Kramers degeneracy can be removed by 
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adding an external orbital-polarizing field, 

h e = -uaL, (36) 

where wa is the cranking frequency for the orbital mo- 
tion. Indeed, in the absence of the spin-orbit coupling 
the spin and the orbital angular momentum may rotate 
with different angular velocities. In the presence of the 
field (l36|) . each level is shifted by —oj\A, i.e., the energy 
splitting of the Kramers doublet becomes 2wa|A|. An 
illustrative example of such situation is displayed in the 
bottom panels of Fig. 3. Here, each level corresponds 
to a one-quasiparticle excitation. We confirmed numeri- 
cally that the result of calculations for N—31 by explic- 
itly blocking the lowest level are here equivalent to those 
with 2FLA carried out above the crossing. It is seen in 
Fig. 3, however, that because of high density of quasipar- 
ticle levels, the self-consistent calculations in 2FLA are 
difficult due to many consecutive crossings that make it 
extremely difficult to keep track of the fixed configura- 
tion. We note that the order of the quasiparticle levels 
in 7V=30 and 7V=31 systems is affected by the variation 
of the mean field due to the crossing. 

In order to illustrate 2FLA in the nuclear case, we car- 
ried out nuclear DFT calculations using the Skyrme en- 
ergy density functional SLy4 [36| in the p-h channel, aug- 
mented by the "mixed-pairing" [37j density-dependent 
delta functional in the p-p channel. The details per- 
taining to the numerical details, ej^ the pairing space 
employed, can be found in Ref. 35]. As a representa- 
tive example, we took the pair of deformed nuclei 166 Er 
(iV=98) and 167 Er (N=99). The pairing strength V =- 
320 MeVfm 3 was slightly enlarged to prevent pairing 
from collapsing in the A^=99 system. The resulting neu- 
tron pairing gaps, A„=1.2MeV and 0.77MeV in 166 Er 
and 167 Er, respectively, are reasonably close to the ex- 
perimental values of 1.02 MeV and 0.62 MeV. 

Figure 4 displays the quasiparticle spectrum for 166 Er 
(top) and 167 Er (bottom). At the value of X s indicated by 
a star symbol, a transition from a zero-quasiparticle vac- 
uum corresponding to A^=98 to a one-quasiparticle vac- 
uum associated with A^=99 takes place. As in the atomic 
case, the mean-field changes abruptly at the crossing 
point. Actually, since the quasiparticle spectrum changes 
when going from 166 Er to 167 Er (both in terms of exci- 
tation energy and ordering of levels), the crossing point 
is shifted towards the lower values of X s . As checked nu- 
merically, the result of calculations for 167 Er by explicitly 
blocking the level "a" are equivalent to those with 2FLA 
carried out above the crossing point. 

V. CONCLUSIONS 

In this study, we analyzed various approaches to po- 
larized Fermi systems within the DFT. The main conclu- 
sions can be summarized as follows: 

• By analogy with rotating nuclei, we showed that 
introducing two chemical potentials for different 
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X s (MeV) 

FIG. 4: (Color online) One-quasiparticle levels of both signa- 
tures (a=l/2: solid line; a=-l/2: dotted line) as functions of 
A s for 166 Er (top) and 167 Er (bottom). The levels occupied 
in the vacuum configuration are drawn by thick lines. The 
calculations were carried out with a SLy4 Skyrme functional 
and mixed pairing. See text for details. 

superfluid components is equivalent to applying a 
one-body, time-odd field. This field can be used to 
change the particle-number parity of the underly- 
ing quasiparticle vacuum. 

• Since the external one-body field commutes with 
the HFB Hamiltonian, 2FLA is equivalent to non- 
collective cranking, a technique that is often used 
to select a vacuum configuration of interest. 

• For systems, in which no additional degeneracy is 
present beyond the Kramers degeneracy, the 2FLA 
is equivalent to one-dimensional, non-collective ro- 
tational cranking. Different choices of the angu- 
lar momentum quantization axis give rise to differ- 
ent blocking procedures and different polarization 
schemes. By increasing the asymmetry A s , one is 
alternating between even and odd systems while 
gradually increasing the signature polarization. In 
the absence of the spin-orbit coupling, the signature 
quantum number can be replaced by spin projec- 
tion. 

• The generalization of 2FLA to the case of an arbi- 
trary blocked state is not obvious, although one can 
introduce polarizing fields that would single out the 
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state of interest. For instance, one can introduce 
different cranking frequencies for spin and orbital 
motion, and this removes the orbital degeneracy of 
quasiparticle states. However, the standard way of 
blocking a given level in a fixed Hamiltonian sub- 
matrix should work as well, and it is easy to gener- 
alize to an arbitrary state, even if the majority of 
self-consistent quantum numbers are gone. 

• For systems, in which an additional fc-fold degener- 
acy (fc-even) is present beyond the Kramers degen- 
eracy, the 2FLA is unable to produce odd-particle- 
number states. Moreover, for an fc-fold degeneracy 
with odd k, which is the case, e.g., for spherical sys- 
tems, the 2FLA gives odd-particle-number states 
that correspond to fc-quasiparticle and not to one- 
quasiparticle excitations. 

• The situation encountered in Fig. 2 is similar to the 
level crossing discussed in the context of high-spin 
physics, where the lowest quasiparticle Routhian 
becomes negative as a function of rotational fre- 
quency. In such situations, one needs to preserve a 
number of quasiparticles in each signature block to 
conserve ttn in the HFB vacuum [23. [3l|. 

• Irrespective of the degeneracies present in the sys- 
tem, certain quasiparticle configurations cannot be 
approached through 2FLA. 

In summary, the 2FLA provides a unifying methodol- 
ogy to treat a number of different kinds of condensates, 



including those of odd-particle systems, as ground states 
of some HFB Hamiltonian. However, not all quasipar- 
ticle states are easily accessible this way and problems 
arise if quasiparticle levels show degeneracies beyond the 
Kramers doubling. Moreover, the examples shown in 
Figs. 3 and 4 indicate that the variations of the self- 
consistent mean field associated with the configuration 
change driven by the polarizing field can be severe, in- 
cluding the change of ordering of the lowest quasiparticle 
excitations. This, together with related numerical in- 
stabilities, indicates that traditional methods of blocking 
are likely to be preferred for treating a large space of 
quasiparticle configurations and in spectroscopic-quality 
calculations for well-defined one-quasiparticle states. 
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